Mon. Not. R. Astron. Soc. 000, [T]-?? (2011) 



Printed 23 November 2011 



(MN XAT^ style file v2.2) 



Recurrent dust formation by WR 48a on a 30-year 
timescale"^ 



o 



Peredur M. Williams^f, Karel A. van der Hucht^'^, Francois van Wyk 



4 



^ Fred Marang'^, Patricia A. Whitelock^'^, Patrice Bouchet^ and 
> 
O 



Diah Y. A. Setia Gunawan^ 

^ Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ 
^ Space Research Organization Netherlands, Sorbonnelaan 2, NL-3584 OA Utrecht, The Netherlands 



^Nj ' ^ Astronomical Institute 'Anton Pannekoek', University of Amsterdam, Science Park 904, NL-1098 XH Amsterdam, The Netherlands 
^Nj : South African Astronomical Observatory, P.O. Box 9, 7935 Observatory, South Africa 

^Astronomy, Cosmology and Gravity Centre, Astronomy Department, University of Cape Town, 7700 Rondebosch, South Africa 
' I ^Service dAstrophysigue DSM/IRFU/SAp CEA - Saclay, L'Orme des Merisiers bdtiment 709, F-91191 Gif-sur-Yvette France 
plL^ ^ Atacama Large Millimetre/ submillimetre Array (ALMA), Alonso de Cordova 3107, Vitacura, Santiago 763 0355, Chile 

in ■ 

^ ^, Accepted 2011 November 16. Received 2011 November 16 ; in original form 2011 August 2 

o 

, ABSTRACT 

52 ' We present infrared photometry of the WC8 Wolf-Rayet system WR 48a observed with 

' telescopes at ESO, the SAAO and the A AT between 1982 and 2011 which show a slow 

decline in dust emission from the previously reported outburst in 1978-79 until about 
1997, when significant dust emission was still evident. This was followed by a slow 
^ \ rise, accelerating to reach and overtake the first (1978) photometry, demonstrating 

i that the outburst observed in 1978-79 was not an isolated event, but that they recur 

at intervals of 32+ years. This suggests that WR48a is a long-period dust maker 
and coUiding-wind binary (CWB). The locus of WR48a in the (Jl — L), K colour- 
magnitude diagram implies that the rate of dust formation fell between 1979 and 
about 1997 and then increased steadily until 2011. Superimposed on the long-term 
variation are secondary ('mini') eruptions in (at least) 1990, 1994, 1997, 1999 and 
2004, characteristic of relatively brief episodes of additional dust formation. Spectra 
show evidence for an Oe or Be companion to the WC8 star, supporting the suggestion 
that WR 48a is a binary system and indicating a system luminosity consistent with the 
association of WR48a and the young star clusters Danks 1 and Danks 2. The range 
^ , of dust formation suggests that these stars are in an elliptical orbit having e ^ 0.6. 

■ The size of the orbit implied by the minimum period, together with the WC wind 

velocity and likely mass-loss rate, implies that the post-shock WC wind is adiabatic 
throughout the orbit - at odds with the observed dust formation. A similar conflict is 
observed in the 'pinwheel' dust-maker WR112. 

Key words: stars: Wolf-Rayet — stars: circumstellar matter — stars: binary — 
stars: individual: WR48a — stars: individual: WR112 
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1 INTRODUCTION 

Population I Wolf-Rayet (WR) stars are characterised 
by dense, fast stellar winds giving the stars their broad 
emission-line spectra and carrying heavy mass loss: they rep- 
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resent the last stable stage in the evolution of massive stars. 
One of the earliest results from infrared (IR) astronomy 
was the discovery of 'excess' radiation by heated circum- 
stellar dust from a variety of stars known to be losing mass. 
Amongst these were three WC9 type WR stars, observed 
by Allen, Harvey & Swings (1972). Subsequent studies (e.g. 
Cohen, Barlow & Kuhi 1975) found circumstellar dust emis- 
sion from most WC9 and many WC8 stars. The process was 
soon inverted, and bright IR sources discovered in space- 
or ground-based IR surveys were found to harbour optically 
faint dusty WC8-9 stars (e.g. AFGL 2104 = WR 112, Cohen 
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& Kuhi 1977). Amongst these was WR48a, discovered by 
Danks et al. (1983, hereafter DDWS) in a survey of the H ll 
region G305.5 +0.0. They observed WR48a to be undergo- 
ing an infrared outburst similar to that shown in 1977 by HD 
193793 (WiUiams et al. 1978, Hackwell, Gehrz & Grasdalen 
1979). DDWS tracked the rise of the IR emission to maxi- 
mum in 1979 and its initial fading, and presented a spectrum 
showing He and C emission lines, most conspicuously He I 
10830A and Cm 9710A, from which they classified WR48a 
as a WC9 star. 

In their survey of IR dust emission by late subtype 
WC stars, Williams, van der Hucht & The (1987a, hereafter 
WHT), reported photometry of WR48a in 1982-86 which 
showed continuing fading of its dust emission. They also ex- 
tended the spectroscopy to the yellow (5550-5950 A), and 
argued from the approximately equal intensities in the C in 
5696- A and C iv 5808- A classification lines that WR48a was 
a WC8 star. Further IR photometry of WR 48a was reported 
by Williams et al. (2003, hereafter Paper I), showing that 
the fading continued to about 1997, but that the spectral en- 
ergy distribution (SED) at that time still showed dust emis- 
sion, indicating that dust was still being formed by WR 48a. 
In this respect it differed from 'episodic' dust makers like 
HD 193793 (WR140), which form dust for only a fraction 
of their periods (Williams et al. 1990). The IR photometry 
of WR48a also showed secondary dust-formation episodes, 
'mini eruptions', in 1990 and 1994. These resembled those 
observed by Williams et al. (2001) in another episodic dust 
maker, the WC7-h09 system WR137 (HD 192641). In Pa- 
per I, the 2.4-25-/im SED of WR48a observed wfth the 
short-wavelength spectrometer (SWS) on ISO in a survey 
of WR spectra (van der Hucht et al. 1996) was compared 
with that of WR104 (Ve2-45), the classical dust-making 
WC9 star, and found to be significantly cooler. This differ- 
ence was consistent with the radial dust density distribution 
around WR48a being flatter than the oc law expected 
of a system making dust at a constant rate like WR 104: 
the evolution of the dust density distribution of WR48a is 
considered further in Section 13.41 The /SO-SWS spectrum 
of WR 48a presented in Paper I showed an emission feature 
near 7.7 /um similar to that observed in WR 104 and WR 112 
by Cohen, Tielens & Bregman (1989) and ascribed to a car- 
bonaceous carrier with small aromatic domains, a possible 
precurson to carbon grains; an independent study of the 
SWS spectrum by Chiar, Peelers & Tielens (2002) identi- 
fied and discussed emission features at 6.4/im and 7.9/im. 

The enduring interest of dust formation by WR stars is 
the great difficulty in forming dust in such hostile environ- 
ments: close to the stars, the stellar radiation fields would 
heat any dust to above its evaporation temperature but at 
greater distances, where the radiation field is sufficiently di- 
luted, the density of the WR stellar wind is too low to al- 
low homonuclear dust condensations - as recognised early 
on by Hackwell et. al. (1979) in the case of WR 140 and 
discussed by WHT for the WC8-9 stars making dust persis- 
tently. Zubko (1998) found that carbon grains having drift 
velocities relative to the WC8-9 wind could grow via colli- 
sion with carbon ions, but the question of grain nucleation 
remains open. A chemical kinetic study of the formation 
of molecular species and carbon grain precursors in WC9 
winds by Cherchneff et al. (2000) found that only C2 was 
formed in useful amounts, and even that required signifi- 



cantly (factor > 10^) higher densities than expected in the 
regions of isotropic, smooth WR winds where dust grains 
are observed. Understanding the physicochemical process of 
dust formation by WR stars appears to have stalled, but the 
incidence of observed dust formation by WR stars in differ- 
ent systems provides valuable information on the origin of 
the significantly higher densities demanded. 

Large-scale, high-density structures in a WR wind are 
provided by compression of the stellar wind in strong shocks 
formed where the wind of a WR star collides with that 
of a sufficiently luminous companion in a binary system 
- a colliding wind binary (CWB) - as occurs in WR 140. 
Usov (1991) suggested that very high compression factors 
(10^ - 10*) could be produced in WR140 if the heated 
and compressed wind was able to cool efficiently. The link 
between the dust-formation episodes and binary orbit in 
WR 140 is provided by the periodic increases by factor of 
~ 40 of the pre-shock wind density at the wind-collision 
region (WCR) for a brief time around periastron passage 
(Williams 1999) in its very elliptical orbit. Similarly, Lefevre 
et al. (2005) showed that the episodic dust-maker WR 137 
was a spectroscopic binary having a period (13.05 yr) equal 
to that of the dust-formation episodes (Williams et al. 
2001), with periastron passage about a year before maxi- 
mum dust emission. More recently, Williams, Rauw & van 
der Hucht (2009a) showed that the radial velocities of the 
WC5-I-09 system WR 19 were consistent with its being a 
high-eccentricity (e = 0.8) binary, having periastron passage 
shortly before its brief, periodic (10.1-yr) dust-formation 
episodes. 

The connection between the persistent dust-formation 
by the WC9 stars and their movement in binary orbits comes 
from IR images of the dust, which show rotating 'pinwheel' 
structures around WR104, WR98a (IRAS 17380-3031) and 
WR112, interpreted to be dust formed and ejected in a 
stream to one side of a rotating binary system observed at a 
relatively low orbital inclination (Tuthill, Monnier & Danchi 
1999, Monnier, TuthiU, & Danchi 1999, and Monnier et al. 
2007). The observation of non-thermal radio emission from 
the 'pinwheel' systems (Monnier et al. 2002) provides further 
support for their being CWBs; but their stellar components 
have not yet been resolved, so that it has not been possible 
to relate the positions of the stars and WCRs to that of the 
dust. 

Monnier et al. (2007) found the dust emission by 
WR48a to be extended at 2/xm; and the dust cloud was im- 
aged at 12.3/im by Marchenko & Moffat (2007), showing the 
dust emission to be extended E-W, suggesting a fiattened 
structure, possibly a fragmentary spiral, seen at moderately 
high inclination. 

DDWS and Danks et al. (1984) associated WR48a wfth 
two newly discovered galactic clusters, Danks 1 (C 1310- 
624) and Danks 2 (C 1309-624) at a distance between ^ 
3.3 and 4.6 kpc. Bica et al. (2004) derived colour-magnitude 
diagrams from deeper photometry and derived distances to 
the clusters of 3.6 ± 0.5 and 3.4 ± 0.2 kpc respectively. The 
connection between WR48a and the clusters was reinforced 
by Clark & Porter (2004) in their study of star formation in 
the region. They suggested that star formation might have 
been triggered by WR 48a and adopted a distance of 4 kpc 
for the complex. This gives a greater distance to WR48a 
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than that (1.2 kpc, van der Hucht 2001) imphed by its spec- 
tral type, and these estimates need to be reconciled. 

By all accounts, WR48a is heavily reddened, but the 
exact value is not certain. DDWS derived A,, ~ 9.2 from 
the 7000-10000-A colour in their spectrum. Williams, van 
der Hucht & The (1987b) derived A„ ~ 8.3 from the 5000- 
7000-A slope in their spectrum observed in 1986, and noted 
that this was consistent with the strength of the 6284-A 
diffuse interstellar absorption band (DIB) and the relation 
between these quantities from other reddened stars. From 
the ISO SWS spectrum of WR48a, Schutte et al. (1998) 
found the interstellar 'silicate' feature to have rg.7 = 0.52. 
The ratio of visual extinction to 'silicate' optical depth in the 
spectra of WR stars and Cyg 0B2 No. 12 showed significant 
scatter, Av/rg.r = 16.6±4.6. Applying this ratio to the rg.y 
observed in WR48a gives Av ~ 8.6 ± 2.4. (A lower value, 
Av (6.1), was derived from rg.y in Paper I but this used 
a different and probably inappropriate calibration.) Baume, 
Carraro & Momany (2009) observed UBVIc photometry in 
2006, and derived E{B - V) ~ 2.4. For a standard red- 
dening law, this implies a lower extinction and, from their 
observed {V — Ic) = 3.88, an intrinsic {V — Ic) ~ 0.9, which 
is abnormally red for an early type star. If some of the ex- 
tinction occurs in the carbon dust made by WR48a, the 
extinction law may not be standard. For example, Menzies 
& Feast (1997) found an anomalous reddening law for the 
carbon dust made by the RGB star RY Sgr, so the sightline 
to WR 48a may have a reddening law intermediate between 
the standard interstellar and that for carbon dust. On the 
other hand, the circumstellar component of the reddening 
is unlikely to be very great, given the strength of the 'sil- 
icate' feature and the reddening to the associated clusters 
Av = 8.1 and 9.3 respectively (Bica et al.). For the present, 
therefore, we adopt Av = 9.0 and a standard interstellar 
law for the reddening of WR48a. 

Zhekov, Gagne & Skinner (2011) have identified a 
strong, highly obscured X-ray source with WR48a. They 
pointed out that, if WR48a was physically associated with 
the clusters Danks 1 and Danks 2 at a distance of ~ 
4kpc, it would be the most X-ray luminous WR star in 
the Gala:xy, after Gyg X-3. They also derived a high ex- 
tinction, Nh = 2.29 X 10^^ cm"^, from their NEI shock 
modelling of the X-ray spectrum and noted that this was 
consistent with Av ~ 10 using the standard conversion, 
Nh = 2.22 X lO'^'^Av cm"^; but Av could be 40-50% higher 
using more recent conversions. As one of several possible ori- 
gins for the X-ray emission, they considered colliding stellar 
winds in a wide WR+0 binary. As noted by Zhekov et al., 
and also Mauerhan, van Dyk & Morris (2011), WR48a was 
detected with Chandra as a very strong X-ray source. 

In this paper, we report almost three decades (1982- 
2011) of IR photometry tracking the emission by the dust 
cloud made by WR48a and showing that the 1978-9 dust- 
formation was not an isolated event, but one which appears 
to be recurring after an interval of over 32 yrs. We present 
optical spectra which suggest that the WG8 star in WR 48a 
has an OB companion, apparently an Oe or Be star. We use 
the IR colours and simple modelling of the SEDs to describe 
the evolution of the dust cloud and develop a preliminary 
picture of WR 48a as a colliding- wind binary. 



Table 1. Near-infrared photometry of WR48a with ESO tele- 
scopes: except where marked (* 3.6m, + 2.2m), observations were 
made with the 1-m photometric telescope. 



Date 


J 


H 


K 


L' 


M 


^'-t'-tO^-LO 


7.80 


5.21 


3.01 


0.29 


-0 42 


^'-t'-to^-L 


7.88 


5.25 


3.10 


0.41 


-0.37 




8.27 


5.61 


3.39 


0.52 


-0.31 




8.46 


5.94 


3.71 


0.80 


04 




8.57 


6.23 


4.07 


1.13 


27 




8.68 


6.40 


4.22 


1.18 


0.32 




8.69 


6.38 


4.12 


1.27 


0.25 


9zL479f^« 


8.72 


6.46 


4.35 


1.41 


58 


9zL479fiQ 


8.73 


6.45 


4.32 


1.36 


63 


9A479fin 


8.70 


6.45 


4.33 


1.40 


0.63 


9A479fi1 


8.70 


6.46 


4.36 


1.42 


0.61 


9zL47f;i Q 


8.91 


6.69 


4.57 


1.54 


64 




8.94 


6.69 


4.57 


1.54 


73 


9Ad7'^Af) 


9.17 


6.78 


4.56 


1.59 


0.67 




8.94 


6.75 


4.55 


1.56 


0.66 


^AdV^Ad 


8.93 


6.70 


4.56 


1.58 


69 


9zL47'^fi9 


8.85 


6.73 


4.55 


1.55 


70 


2447^73 


8.87 


6.71 


4.60 


1.58 


0.74 


9447(^70 


8.73 


6.67 


4.63 


1.69 


0.91 


2447671 


8.86 


6.80 


4.70 


1.66 


0.73 


9zL47fi79 


8.86 


6.82 


4.73 


1.71 


78 


2447910 


8.89 


6.83 


4.80 


1.79 


0.89 




8.87 


6.87 


4.86 


1.85 


0.94 


9447Q(S7 


8.93 


6.89 


4.86 


1.87 


0.96 




8.90 


6.87 


4.83 


1.83 


92 


2447970 


8.91 


6.88 


4.85 


1.84 


0.96 


9A47Q71 


8.90 


6.88 


4.84 


1.83 


0.95 


2447Qqi 


8.91 


6.89 


4.87 


1.87 


0.93 


9447QQ9 


8.86 


6.88 


4.86 


1.84 


95 


2447QQS 


8.89 


6.88 


4.87 


1.85 


0.97 


244SS1 fi 


8.89 


6.91 


4.89 


1.94 


1.01: 




8.89 


6.89 


4.89 


1.92 


1.02 




8.91 


6.90 


4.88 


1.93 


1 02 


9zL4S'?zL7* 


8.91 


6.93 


4.93 


1.97 


1 07 


'?AdR'^A^ 


8.94 


6.96 


4.95 


1.97 


1.04 




8.96 


6.96 


4.96 


1.98 


1.09 


244S3'in 


8.95 


6.98 


4.97 


1.97 


1.08 


2448639 


8.91 


6.92 


4.95 


2.03 


1.16 


2448698 


8.96 


7.02 


5.04 


2.08 


1.14 


2448749 


8.98 


7.07 


5.10 


2.11 


1.21 


2448752 


9.00 


7.07 


5.10 


2.12 


1.20 


2448753 


8.92 


7.05 


5.09 


2.12 


1.21 


2448799 


9.16 


7.14 


5.13 


2.13 


1.20 


2449050+ 


9.01 


7.15 


5.24 


2.26 


1.35 


2449052+ 


8.98 


7.14 


5.24 


2.26 


1.35 


2449179 


8.97 


7.09 


5.20 


2.29 




2449180 


8.96 


7.10 


5.21 


2.30 


1.37 


2449181 


8.94 


7.11 


5.22 


2.31 


1.38 


2449182 


8.92 


7.11 


5.23 


2.32 


1.39 


2449183 


8.96 


7.11 


5.22 


2.31 


1.36 


2450198+ 




7.11 


5.37 


2.49 


1.61 


2450890+ 


9.07 


7.12 


5.20 


2.54 


1.78 



2 OBSERVATIONS 

2.1 Infrared photometry 

Photometry of WR48a was obtained using aperture pho- 
tometers on the 1-m, 2.2-m and 3.6-m telescopes of the Euro- 
pean Southern Observatory (ESO) at La Silla, Ghile, (Tables 
[T]and[2}, the Mk II photometer on 0.75-m telescope of the 
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Table 3. SAAO photometry of WR48a with the Mk II IR photometer on the 0.75-m telesope using the 36"(*or 24") aperture. 



Date 


J 


H 


K 




L 


2445417 


7.78 


5, 


.19 


3, 


,05 


0, 


,65 


2447547 


8.38 


6, 


,50 


4, 


,50 


1, 


,98 


2447553 


8.34 


6, 


.49 


4, 


,57 


1, 


,98 


2447607 


8.34 


6, 


.55 


4, 


,57 


2, 


,00 


2447634 


8.42 


6, 


.64 


4, 


,63 


2, 


,11 


2447669 


8.39 


6, 


,59 


4, 


,62 


2, 


,08 


2447687 


8.43 


6, 


,63 


4, 


,63 


2, 


,08 


2447737 


8.42 


6, 


,65 


4, 


,67 


2, 


,17 


2448028 


8.36 


6, 


,66 


4, 


,77 


2, 


,24 


2448108 


8.31 


6, 


,46 


4, 


,61 


2, 


,30 


2448285 


8.40 


6, 


,87 


4, 


,78 


2, 


,29 


2448289 


8.42 


6, 


,70 


4, 


,78 


2, 


,31 


2448327 


8.41 


6, 


,73 


4, 


,82 


2, 


,34 


2448365 


8.43 


6, 


,76 


4, 


,87 


2, 


,35 


2448378 


8.38 


6, 


,73 


4, 


,85 


2, 


,41 


2448385 


8.41 


6, 


,76 


4, 


,87 


2, 


,35 


2448427 


8.42 


6, 


,74 


4, 


,88 


2, 


,42 


2448437 


8.41 


6, 


,74 


4, 


,89 


2, 


,41 


2448451 


8.40 


6, 


,73 


4, 


,86 


2, 


,39 


2448457 


8.41 


6, 


,71 


4, 


,85 


2, 


,39 


2448462 


8.40 


6, 


,72 


4, 


,85 


2, 


,38 


2448500 


8.33 


6, 


,62 


4, 


,78 


2, 


,37 


2448701* 


8.47 


6, 


,80 


4, 


,96 


2, 


,48 


2448734* 


8.54 


6, 


,83 


4, 


,98 


2, 


,51 


2448760* 


8.44 


6, 


,83 


5, 


,00 


2, 


,51 


2448783 


8.40 


6, 


,82 


5, 


,00 


2, 


,52 


2448792 


8.46 


6, 


,85 


5, 


,01 


2, 


,54 


2448817 


8.45 


6, 


,86 


5, 


,03 


2, 


,56 


2448851 


8.51 


6, 


,90 


5, 


,07 


2, 


,63 


2449001 


8.48 


6, 


,94 


5, 


11 


2, 


,60 


2449106 


8.52 


6, 


,91 


5, 


14 


2, 


,68 


2449146 


8.42 


6, 


,87 


5, 


12 


2, 


,69 


2449221 


8.44 


6, 


,90 


5, 


14 


2, 


,75 



Date 


J 


H 


K 


L 


2449382 


8.56 


6.97 


5.23 


2.79 


2449406 


8.46 


6.86 


5.14 


2.77 


2449501 


8.35 


6.62 


4.91 


2.68 


2449519 


8.28 


6.55 


4.85 


2.67 


2449583 


8.32 


6.58 


4.84 


2.62 


2449590 


8.34 


6.59 


4.84 


2.59 


2449735 


8.43 


6.75 


4.96 


2.61 


2449820 


8.46 


6.79 


4.98 


2.68 


2449887 


8.46 


6.85 


5.04 


2.68 


2450091 


8.47 


6.91 


5.13 


2.76 


2450127 


8.44 


6.89 


5.13 


2.77 


2450178 


8.45 


6.95 


5.21 


2.75 


2450206 


8.42 


6.95 


5.20 


2.83 


2450237 


8.47 


6.95 


5.22 


2.86 


2450259 


8.49 


6.95 


5.23 


2.90 


2450320 


8.48 


6.99 


5.28 


2.90 


2450484 


8.46 


7.02 


5.32 


2.90 


2450591 


8.46 


7.07 


5.44 


3.03 


2450618 


8.52 


7.06 


5.43 


3.07 


2450798 


8.41 


6.82 


5.12 


2.89 


2450980 


8.44 


6.92 


5.23 


2.98 


2451234 


8.40 


6.83 


5.14 


2.92 


2451298 


8.46 


6.93 


5.23 


2.94 


2451358 


8.45 


6.93 


5.21 


2.99 


2451576 


8.34 


6.68 


5.02 


2.89 


2451614 


8.39 


6.71 


4.99 


2.86 


2451677 


8.42 


6.75 


5.02 


2.86 


2451685 


8.45 


6.77 


5.02 


2.87 


2451715 


8.42 


6.79 


5.05 


2.84 


2451743 


8.43 


6.82 


5.06 


2.86 


2451785 


8.44 


6.84 


5.08 


2.88 


2451967 


8.42 


6.86 


5.14 


2.93 



Date J H K L 



2451981 


8.42 


6 


.86 


5, 


14 


2, 


,89 


2452070 


8.40 


6 


.78 


5, 


,06 


2, 


,90 


2452093 


8.41 


6 


.79 


5, 


,06 


2, 


,87 


2452116 


8.45 


6 


.80 


5, 


,07 


2, 


,88 


2452290 


8.52 


6 


.85 


5, 


10 


2, 


,90 


2452322 


8.41 


6 


.81 


5, 


,08 


2, 


,91 


2452351 


8.41 


6 


.84 


5, 


12 


2, 


,90 


2452383 


8.42 


6 


.80 


5, 


11 


2, 


,94 


2452426 


8.36 


6 


.75 


5, 


,07 


2, 


,93 


2452451 


8.38 


6 


.74 


5, 


,05 


2, 


,91 


2452694 


8.44 


6 


.87 


5, 


18 


2, 


,97 


2452755 


8.36 


6 


.72 


5, 


,03 


2, 


,92 


2452798 


8.55 


6 


.82 


5, 


,05 


2, 


,94 


2452868 


8.38 


6 


.68 


4, 


,97 


2, 


,86 


2453065 


8.41 


6 


.72 


4, 


,92 


2, 


,76 


2453118 


8.36 


6 


.66 


4, 


,90 


2, 


,74 


2453177 


8.22 


6 


.42 


4, 


70 


2, 


,61 


2453194 


8.24 


6 


.39 


4, 


,66 


2, 


,60 


2453408 


8.26 


6 


.37 


4, 


,54 


2, 


,43 


2453442 


8.31 


6 


.42 


4, 


,55 


2, 


,42 


2453507 


8.30 


6 


.43 


4, 


,60 


2, 


,43 


2453562 


8.28 


6 


.38 


4, 


,55 


2, 


,45 


2453906 


8.19 


6 


.05 


4, 


19 


2, 


,14 


2453935 


8.13 


6 


.06 


4, 


,20 


2, 


,13 


2454184 


7.99 


5 


.77 


3, 


,88 


1, 


,84 


2454234 


7.93 


5 


.70 


3, 


,81 


1, 


,75 


2455011 


7.52 


5 


.04 


3, 


16 


1, 


,08 


2455353 


7.56 


4 


.90 


2, 


,99 


0, 


,99 


2455621 


7.21 


4 


.71 


2, 


,86 


0, 


,82 


2455642 


7.25 


4 


.75 


2, 


,92 


0, 


,99 


2455731 


7.26 


4 


.73 


2, 


,84 


0, 


,77 


2455733 


7.24 


4 


.72 


2, 


,85 


0, 


,77 



Table 2. Long-wavelength photometry of WR48a with ESO tele- 
scopes 



Date 


L' 


M 


Nl 


N2 


N3 


QO 


Tel. 


2445150 


0.08 


-0.67 


-1.6 


-1.4 


-1.9 


-1.7 


3.6m 


2445425 


0.29 


-0.46 


-1.4 


-1.3 


-1.6 


-1.5 


3.6m 


2446194 


0.86 


-0.03 


-1.3 


-1.1 


-1.7 




3.6m 


2447519 


1.54 


0.73 


-0.9 


-0.5 


-1.4 




1.0m 


2447987 


1.81 


0.80 


-0.8 


-0.6 


-1.2 


-1.2 


2.2m 


2447988 


1.88 


0.91 


-0.7 


-0.5 


-1.3 


-1.3 


2.2m 


2448346 


1.95 


1.10 


-0.6 


-0.4 


-1.1 


-1.2 


3.6m 


2448347 


1.96 


1.09 


-0.6 


-0.3 


-1.1 


-1.1 


3.6m 



South African Astronomical Observatory (SAAO), Suther- 
land, South Africa (Table [3]), and the Infrared Photome- 
ter/Spectrometer (IRPS) on the Anglo-Australian Telescope 
(AAT) at Siding Spring, Austraha (Table gj. The final two 
observations in Table[T]were obtained with the array camera 
IRAC-1 on the 2.2-m telescope at La Silla. 

During our campaign, the near-infrared (NIR) pho- 
tometric systems on the three ESO telescopes were the 
same, and as described by Bouchet, Manfroid & Schmider 
(1991) and Bersanelh, Bouchet & Falomo (1991). The final 
('Postscript') photometry by DDWS was on this system, 
whereas that in the body of their paper was on the earlier 



ESO system (Wamsteker 1981), which differed principally in 
using a filter having Acff — 3.5-3.6 fim instead of 3.8 fim for 
L. The 3.8-^m filter is often called 'L'' to distinguish it from 
the shorter-wavelength L, a convention we follow here. The 
SAAO photometry (Carter 1990) uses a shorter-wavelength 
(Acff ~ 3.5Atm) L filter. The NIR system at the AAT (Allen 
& Cragg 1983) also used a 3.8-pim filter for L' , but their J 
filter had a slightly shorter Aofj than the ESO J filter and 
closer to that of the J filter in the SAAO Mk II photometer 
(see Glass 1993). 

Observations of WR 48a can be affected by the nearby 
11th magnitude star CPD -62° 3058. From HST WFPC2 im- 
ages, Wallace et al. (2003) measured a separation 9'.'2 and 
position angle 194? 2, placing it within or near the edges of 
some of the apertures used for the photometry of WR 48a. 
To determine the effect of this contribution on the photom- 
etry, we observed CPD -62° 3058 with IRIS on the AAT 
through a 5" aperture and give the results, together with 
magnitudes from the DENIS (Epchtein et al. 1999) and 
2MASS (Skrutskie et al. 2006) surveys, in Table [3 Despite 
the scatter in the data, we assume that the CPD star is 
not significantly variable in the NIR. In the J band, the 
CPD star is 0.2 mag. fainter than WR48a at its minimum 
(2.4 mag. at maximum) and so makes a significant contri- 
bution to the J magnitudes observed with the SAAO 0.75- 
m telescope, which used 36" or, occasionally, 24" apertures 
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Table 4. Infrared photometry of WR48a with the IRPS on the 
AAT 



Date 


J 


H 


K 


L' 


M 


2446309 


8.61 


5.94 


3.73 


1.00 


-0.11 


2446595 


8.88 


6.17 


4.00 


1.11 




2446716 


8.98 


6.28 


4.07 


1.18 




2446888 


8.76 


6.40 


4.26 


1.31 




2447614 


8.97 


6.77 


4.70 


1.71 





Table 5. Infrared photometry of CPD -62° 3058 with the IRPS 
or from surveys. 



Date 


J 


H 


Ks 


K 


L' 


source 


2446888 


9.48 


9.20 




8.99 




AAT 


2447614 


9.39 


9.13 




8.97 


8.53 


AAT 


2450501 


9.33 




8.86 






DENIS 


2450973 


9.28 










DENIS 


2451595 


9.30 


9.13 


8.90 






2MASS 



(see Table [3]) . This contribution is much smaller at longer 
wavelengths, with differences (WR48a minus CPD star) in 
the ranges AH = 2 — 5 mag., AK = 3.5 — 6.2 mag., and 
AL ~ 5 — 8 mag. corresponding to minima and maxima in 
WR48a. The AAT photometry used apertures of 5" or 6" 
and that with the ESO 3.6-m telescope used a 7'.'5 aper- 
ture, so the CPD star should have been well excluded. De- 
pending on the seeing and centring, there might have been 
some contamination of observations with the ESO 1-m and 
2.2-m photometers, which generally used 15" and 12" aper- 
tures. We searched for this by comparing the dispersion in 
the J magnitudes when WR 48a was observed on successive 
nights in an observing run (Table [1} with the dispersion in 
K, where any contamination by the CPD star would be at 
least a factor of 20 smaller, but found no significantly greater 
dispersion in J. 

2.2 Homogenizing the NIR photometry 

In order to determine the evolution of the dust emission 
by WR48a from the IR photometry, it is necessary to con- 
sider the relations between the different systems. The zero 
points of the systems are fairly similar (Table |6l compiled 
from Cohen et al. 1999) but the redness of WR48a com- 
pared with the photometric standards is more of a problem. 
Bouchet et al. (1991), Carter (1990), and Allen & Cragg 
(1983) all derive colour transforms between their and other 

Table 6. Zero points for the photometric systems employed 
and our ad hoc adjustments appUed to the L/L' photometry of 
WR48a, the bands being identified by their wavelengths. 



System 


H 


K 


3.5 


3.6 


3.8 


AAT 


-0.02 


-0.01 






-0.03 


pre-1982 ESO 


-0.03 


-0.01 




0.04 




post-1982 ESO 


-0.06 


-0.03 






-0.01 


SAAO 


0.00 


0.00 


0.00 






adjustments 






0.00 


0.18 


0.38 



NIR photometric systems, none of these extend to objects 
as red as WR 48a: most comparisons are based on stars hav- 
ing {J — K) < 1.1, whereas WR48a generally has (J — K) 
~ 4.0. We therefore estimated ad hoc transformations for 
WR48a, concentrating on the various L and L' systems, 
as the filters for these have very different wavelengths. We 
relate the AAT and two sets of ESO photometry to the cor- 
rected SAAO data as the last forms the largest homogeneous 
dataset. 

To relate the post-1982 ESO photometry of WR48a to 
the corrected SAAO photometry, we compared contempora- 
neous observations, i.e., those made within about 0.02 yr of 
each other, sometimes using successive observations at one 
observatory to bracket those made at the other and avoid- 
ing times when WR48a appeared to be changing rapidly. 
This gave eight useful comparisons between 1983 and 1993, 
yielding differences in the sense SAAO minus ESO: AJ = 
-0.09 ± 0.04, AH = 0.07 ± 0.03, AK = 0.04 ± 0.01 and 
L ~ L' = 0.38 ± 0.06, while the ESO (J - K) colour of 
WR48a was in the range 4.8-3.9. The difference L — L' is 
much greater than that found by Bouchet et al. (1991) for 
stars having (J — A")<1.0, but not surprising given the dif- 
ference between the two filters and the extreme redness of 
WR48a. Attempts to derive a colour equation, such as a lin- 
ear fit of the differences to (J — K), did not yield statistically 
significant coefficients, so we restricted ourselves to applying 
a constant difference, L — L' = 0.38, to our ESO photometry 
(and to the final, 1982, observation by DDWS) of WR48a 
to bring the data into line with the SAAO photometry for 
the production of homogeneous light and colour curves. 

For the earlier ESO photometry of WR48a in DDWS, 
we do not have contemporaneous SAAO observations for di- 
rect comparison. Instead, we used the photometry of WR 
stars by Pitault et al. (1983), which was also on the older 
ESO system and included two dust-emitting WC9 stars hav- 
ing (J — K) comparable to that of WR48a as an interme- 
diate step. Comparison of their photometry of WR 104 and 
WR 106 with observations of these stars by WHT, which 
was on the post-1982 ESO system, suggests an offset L-L' — 
0.20. This is about half the difference between L'and the 
SAAO L, as one might expect from the fact that the efec- 
tive wavelength of the older ESO L filter lies between those 
of L'and the SAAO L. To convert DDWS L magnitudes of 
WR48a to the SAAO L scale, we therefore added 0.18 mag., 
and made no adjustments to H or K. 

Four of the AAT observations of WR48a are contem- 
poraneous with observations at ESO or SAAO, from which 
we infer an offset L-L'= 0.34, comparable to that for the 
ESO photometry, as one might expect from the similarity of 
the filter wavelengths. Inspection of the AAT H and K data 
suggests that there is no significant shift between them and 
the ESO data. 

We applied these differences to the AAT and the two 
sets of ESO photometry and combined them with the cor- 
rected SAAO photometry to produce a homogeneous pho- 
tometry set in HKL. The light curves from these data and 
the mid-IR photometry in Table [2] are discussed below, and 
the trajectory in the near-IR colour-magnitude diagram is 
discussed in Section [3.21 The adjusted L/L' magnitudes are 
not used for production of the SEDs in Sections 13.31 and 13.41 
where the calibrated fluxes are entered at their respective 
wavelengths. 
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Figure 1. Optical spectra of WR48a. Top: 18-A-resolution spec- 
trum observed in 1993. The 10830-A line is truncated to allow 
the weaker lines to be more easily seen. Middle: 7.4-A-resolution 
spectrum observed with the blue arm on the same time. Bottom: 
blue spectrum of WR48a observed in 2001 with EMMI on the 
NTT, blocked to 3-A bins. 



2.3 Optical spectroscopy 

New spectra were observed for us with the RGO spectro- 
graph on the Anglo Australian Telescope on 1993 June 21, 
when the NIR emission had faded about 2.5 mag. from max- 
imum (cf. Fig[3|. In the blue, the resolution was 7.4A while 
that in the red, observed simultaneously, was ISA. A spec- 



trum of the metal-poor star HD 184711 observed on the 
same night was used to help with removal of telluric ab- 
sorption features in the red. Rectified spectra are shown in 
the first two panels of Fig. [1] The correction for the strong 
telluric absorption in the A band (7630A) was incomplete 
because HD 184711 was observed at a rather lower airmass 
than WR48a. 

The 'red' spectrum resembles that in DDWS, but the 
emission lines appear more strongly. We measured equiva- 
lent widths (EWs) of the strong He I 10830-A and Cm 9715- 
A lines to be 280A and IISA respectively, compared with 
EWs of - 63A and 46A estimated from fig. 3 of DDWS. 
Taken at face value, these differences suggest that, at the 
time (1980 January) DDWS's spectrum was observed, these 
emission lines were diluted by additional contimuum emis- 
sion of factors 1.6 - 3.4 above the stellar continuum be- 
tween 9715A and 10830A possibly by the newly formed 
dust, which is considered below. 

The 10830-A emission line falls close to the blue edge 
of the J filter and line emission can contribute to the J 
magnitudes. For example, the transmission of the J filter in 
the SAAO Mk II IR photometer (Glass 1993) is about 50% 
maximum at this wavelength while that of the post-1982 
ESO J filter (Bersanelli et al.) is nearer 20%. The other 
optical components and the atmosphere will reduce these 
percentages. The corresponding contributions of the 10830- 
A line to J magnitudes observed with the Mk II photometer 
at the times of the 1980 and 1993 spectra would then be 
0.01 and 0.05 mag. for a fiat continuum across the J-filter 
passband, and somewhat less for a rising dust continuum, 
while the contributions to the ESO J magnitudes would be 
about half these values. We did not adjust J photometry 
en masse but did adjust the SAAO J magnitude used for 
the SED near minimum (Section 13. 3p for the 10830-A line 
contribution. 

We were not able to measure a P Cygni absorption com- 
ponent of the 10830-A line, but its FWZI of 90 A allows us 
to estimate the wind velocity to be Vod — 1200±170 km s^^. 
This is lower than the average (1700 km s-2) for WC8 stars 
(van der Hucht 2001, table 27) and equal to that for WC9 
stars. This line needs to be re-observed at higher resolution 
and higher S/N ratio to to refine our measurement. (The 
broad C iii feature at 9715A is a multiplet with components 
ranging 9704-9720A.) 

We measured the 6284-A diffuse interstellar band to 
have an EW of 3.1 A, comparable to the value (2. 8 A) in the 
1986 spectrum, observed when the dust emission was greater 
(Williams, van der Hucht & The 1987b) , so it is unlikely 
that this feature has a circumstellar component. 

After the 10830-A line, the strongest emission line is 
that at 6570A (EW 73 A). This is much stronger than the 
C III 9715-A feature, whereas it is weaker than the C ill lines 
in the spectra of the WC8-I-09 systems WRll (7 Vel) and 
WR70 presented by Vreux, Dennefeld & Andrillat (1983), 
and significantly weaker in the single WC8 star in their sam- 
ple, WR53 (since reclassified WC9 by Crowther, de Marco 
& Barlow 1998). The identification of the 6570-A feature 
is difiicult. DDWS identified it with A6560 Hen (Hu/3), but 
we consider any contribution of He 11 - which we cannot rule 
out given the resolution of the spectrum - to be small, be- 
cause of the weakness or absence of the Hua line at 10124A. 
There may be a small Hua contribution to the C lll-iv fea- 
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ture at lOlSOA, but the Hu/3 line would be even smaller. 
Another contributor to the 6570A feature is multiplet 2 of 
Cii. The spectra of WRll, WR53 and WR70 (Vreux et 
al.) suggest that the strength of this should be compara- 
ble (within a factor of 2) to that of multiplet 3 of C ll at 
7236A; but the 6570-A line is a factor of about 16x stronger 
in our spectrum of WR48a, so Cll cannot be the dominant 
contributor. Another possible contributor to the 6570A fea- 
ture in low-resolution spectra of WR48a is the Civ line at 
6592A. This line accounts for the rising ratio of the 6570A 
and 7236A C ll lines in the spectra of earlier subtype WC 
stars (cf. Vreux et al.), but nowhere is the ratio as high as it 
is observed in WR48a - nor are the other Civ lines in our 
spectrum found to be anomalously strong. No other plausi- 
ble WC8 lines seem capable of explaining the strength of the 
6570 A feature. We suggest it comes from Ha, either from a 
nebula very close to the star (not resolved in the long-slit 
spectrum) or from an Oe or Be stellar companion to the WC 
star. 

Owing to the heavy reddening of WR 48a, the blue spec- 
trum is very noisy at short wavelengths, and only that long- 
ward of 4200A is considered useful. The only convincing fea- 
ture in the blue is the blend of C ill 4649A and He ll 4686A. 
This is very weak, having an equivalent width (EW) of 36A, 
compared with the EW of 410A of the same feature mea- 
sured from the low-resolution spectrum of the single WC8 
star WR135 observed by Torres & Massey (1987). There 
is a possible feature at the position of the 4325-A C ill line 
and a stronger, narrow line at 4868 A which might be H/3 
emission from the proposed companion. 

The spectrum in the yellow is similar to that observed 
in 1986 (WHT). Compared with the WC8 spectra in Torres 
& Massey, the Hel 5896-A line in WR48a appears abnor- 
mally strong relative to the C ill and C iv lines. This is not 
a question of the classification of WR 48a: the WC9 spectra 
in Torres & Massey have the He I comparable in strength to 
the Civ line but much weaker than Cm. The carbon lines 
are much weaker: the EWs of the 5696-A and 5813-A fea- 
tures in WR48a are 49 A and 41 A repectively, one quarter 
their values in WR135 (197 A and 185 A respectively). 

We take relative weakness of the C ill and C IV lines as 
evidence for their dilution, in this wavelength region, by the 
proposed stellar companion. Taking the ratios of 5696-A and 
5813-A lines in WR48a and the WC8 standard WR 135, the 
companion would be ~ 3x (1.2 mag.) brighter than the WC8 
star in V and the total luminosity ~ 4x (1.5 mag.) brighter. 
If the WC8 star had My = -3.7 (cf. M„ = -3.7 ± 0.5 
for WC8 stars in clusters and associations, van der Hucht 
2001), the companion would have My — —4.9, compara- 
ble to that of an 06V star (Martins, Schaerer & Hillier 
2005), and the system would have My — —5.2. At 4 pc, 
this would be Vo = 7.8 which, reddening by Av = 9, would 
give V = 16.8, close to the V — 17.1 observed by Baume 
et al. and potentially resolving the discrepancy between the 
spectroscopic distance to WR48a and that of the clusters. 
Fitting the observed V magnitude and retaining the adopted 
distance and reddening would give 05V for the companion 
and Mv ~ —4.0 for the WC8 star, still within the range for 
its spectral type. The total luminosity of the WR48a system 
would be ~ 4 X W^Lq, which matches that found by Clark 
& Porter (2004). 

In an attempt to classify the companion, we re-observed 
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Figure 2. Two-micron spectrum of WR48a observed with IRIS 
and formed from the 5th and 6th order HK echelle spectra. The 
spectrum of WR 48a at this wavelength is dominated by the dust 
emission, most of which is suppressed in the figure to show up the 
Unes. 

WR48a in the blue using the ESO Multi-mode Instrument 
(EMMI) on the 3.6-m New Technology Telescope (NTT) at 
La Silla on 2001 June 17-19 to look for absorption lines. The 
spectrum in the bottom panel of Fig.[l] is blocked down to 
3- A bins to improve the S/N ratio. The Cm 4649- A triplet 
is present about 50% stronger than in the low-resolution 
spectrum but we do not confirm the C ill 4325-A line, nor 
do we see any evidence for emission at the position of C iv 
444lA usually seen more strongly in WC8 and still present in 
WC9 spectra. Evidently, the continuum from the companion 
is dominant in this wavelength region and it is disappointing 
not to observe absorption lines from it, apart from a possible 
weak line at the position of B7, which may owe its weakness 
to partial infill by emission. 

The formal classification of the WC star is clearly WC8 
from the ratio of the C ill 5696-A and C iv 5813-A features 
but it appears to be atypical in terms of Hel strength and 
wind velocity. Without the absorption lines, the companion 
cannot be classified, but estimates of its luminosity adopting 
a distance of 4 kpc and Av = 9 suggests an 05Ve star and 
gives Mv = -4.0 for the WC8 star. 



2.4 Infrared spectroscopy 

The infrared spectrum of WR 48a is dominated by dust emis- 
sion, but it is still possible to see the line spectrum on top 
of this in the near IR. Low-resolution NIR spectra were 
observed for us using the InfraRed Imaging Spectrograph 
(IRIS) on the AAT on 1995 August 15. The HK echelle 
gave a spectral resolution of i? = 400. The spectrum in the 2- 
micron region is shown in Fig.|2] Comparison with the spec- 
tra of WC8 and WC9 stars discussed by Eenens, Williams 
& Wade (1991) and Figer, McLean & Najarro (1997), par- 
ticularly the relative strengths of the Civ triplet at 2.07- 
08 ^m and Clll/Hel blend at 2.11 nm shows that WR48a 
is a closer match to the WC9 spectra than the WC8 star 
WR 135. Stronger than either of these features in WR 48a is 
the Hel 2.06-/im line, also strong in the WC9 stars WR121 
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and WR119 (cf. Figer et al.) but not WR88 (cf. Eenens et 
al.). Of particular interest is the strength of the 2.16-/im He I 
feature, which is stronger then the other lines in any of the 
WC8 or WC9 spectra in Eenens et al. and Figer et al. We 
suggest one possible reason for this is a contribution from 
Br 7 from the companion to the WR star. 



3 EVOLUTION OF THE DUST EMISSION 
FROM WR48A 

3.1 Infrared light curves 

The light curves m H , K and L from the homogenized NIR 
data are shown in Fig.jS] We see that the flux faded un- 
til about 1997, after which it grew very slowly until about 
2003, and then began rising again. In K it reached the level 
of DDWS's first (1978) observation in 2010 and continued 
brightening, albeit more slowly. 

Superimposed on the long-term variation are secondary 
episodes in 1990, 1994-5, 1999-2000, 2004 and perhaps 
1997-8 and 2010-11. The best defined is that in 1994-95. 
From 1994.1, the flux rose to maximum in 1994.4 (J and 
H) to 1994.6 {K and L) and then gradually faded through 
1995 and 1996. The rates of fading are strongly wavelength 
dependent (~ 0.38 mag.y"^ in J, 0.34 mag.y"^ in H, 0.22 
mag.y"^ in K and 0.15 mag.y"'^ in L) and the evolution of 
the IR colour (Section 13. 2p is consistent with the formation 
and cooling of extra dust. Within the coverage of the obser- 
vations, the other episodes have similar shapes: a steep rise 
for about 0.4 yr followed by a slower decline. Simple mod- 
elling of the SED of the epsiode is discussed in Section 13.41 
To determine the timing of other episodes, we shifted the H 
and K light curves defining the 1994 mini outburst to those 
near the other events and used the best fits to measure the 
shifts to the other events, deriving -3.95 yr, -|-3.3 yr, -1-5.4 
yr and -1-10.2 yr, giving intervals between the mini epsiodes 
of 3.95 yr, 3.3 yr, 2.1 yr and 4.8 yr. There is no suggestion 
of periodicity in these data. The J, H and L photometry of 
DDWS suggests that there was another in 1980, but this is 
not evident in their K photometry. It is possible that there 
were mini episodes between then and 1990, but the data are 
too sparse to be certain. As discussed in Section 13.21 the 
trajectory in the NIR colour-magnitude diagram points to 
more secondary outbursts, including a possible one in 2010- 
11. These irregularities suggest that the IR light curves will 
not repeat exactly and will make it hard to determine a 
photometric period without several more years observations 
but the present data suggests a period a little longer than 
32 years and we will discuss WR 48a in the framework of its 
being a long-period CWB. 

We tracked the fading of the dust emission at longer 
wavelengths and show light curves in A'^3 (12.8/^m), A'^1 
(8.4/im), M and L, together with from DDWS in Fig.g] 
As can be seen from Table [2l the N2 and Q magnitudes fol- 
low the same trend (A^2 is always fainter than A^l and 7V3 
owing to the coverage by this band of the strong silicate ab- 
sorption feature) but are omitted from Fig. |4] for clarity. Be- 
tween 1983 and 1991, the rates of fading are strongly wave- 
length dependent: 0.07 mag yr~^ in NZ, 0.08 mag yr~^ in 
A^^l, 0.16 mag yr~^ in M and 0.17 mag yr~^ in L, showing 
continued cooling of the dust. Between the 1979 maximum 
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Figure 4. Mid-IR light curves of WR48a during 1979 maximum 
and subsequent fading. 

and our first mid-IR photometry the fading appears to be 
slower, and at 20 microns our QQ = — 1.7±0.2 is compara- 
ble to O = —1.652 ± 0.2 observed by DDWS, reminiscent of 
the behaviour of WR 140, where the long-wavelength fluxes 
brightened a year after maximum (Williams et al. 2009b). 

3.2 Synopsis of dust evolution from the infrared 
colour-magnitude diagram 

There are three principal processes determining the evolu- 
tion of the dust made by a system like WR48a and its in- 
frared emission (cf. consideration of the dust made by the 
archetypal dust-making CWB WR 140 by Williams et al. 
2009b). These are the nucleation of new grains within the 
compressed stellar winds, the growth of the grains by ac- 
cretion of carbon ions as they are accelerated through the 
wind by radiation pressure, and the cooling of the grains 
as they move away from the stars and the stellar radiation 
heating them is progressively diluted. The nucleation pro- 
cess is not understood, observations indicate that the newest 
grains have temperatures in the region of 1200 K (e.g. WHT, 
Williams et al. 2009b). The growth of the grains can also 
cause them to cool as they become relatively more efficient 
coolants than small grains. Their velocity from the stars is 
initially that inherited from the compressed wind, and is in- 
creased by radiation pressure on the grains until they reach 
a drift velocity relative to the wind when the radiative ac- 
celeration is balanced by supersonic drag. 

Compared with the case of WR 140, we know less about 
the WR48a system, e.g. velocity of the dust and its dis- 
tance from the stars and stellar luminosity, and modelling is 
premature, but we can at least qualitatively identify these 
processes from the SEDs of WR48a and its IR colour and 
magnitude. 
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Figure 3. Synoptic H, K and L light curves of WR48a from SAAO data corrected for contamination from CPD -62° 3058 and the ESO 
and AAT data adjusted as described in the text. Given the long time-scales, we use civil dates for discussion of the dust emission. The 
L magnitudes from different sources are plotted with different symbols and show concordance where they overlap. 



We formed a colour-magnitude diagram from our cor- 
rected and adjusted photometry to track the evolution of the 
dust. For the amount of hot dust, we chose K magnitudes, 
and straddled this with [H — L) as the best indicator of 
the temperature of the newly formed hot dust. Sometimes 
WR48a was observed several times in a short interval, so 
we have gridded the data to avoid confusion in the diagram 
and also increase the accuracy of the points where they come 
from averages of two or more observations. From experimen- 
tation, we found a grid size of 0.2 yr gave a good balance 
preserving information and providing clarity. The CMD is 
shown in Fig.O 

Initially, in 1978-1979.2 (DDWS's observations between 
JD 2443739 and 2443954), we observe the flux brighten- 
ing while the colour did not change. We attribute this to 
the formation of new dust at a higher rate than before. 
The constant colour implies a constant temperature, that 
at which the dust was condensing, and its persistence shows 
that the SED was dominated by emission from the hottest, 
newest dust, - i.e. that the nucleation rate exceeded the 
rate at which the dust was carried away in the wind. Dur- 
ing 1979, the dust emission cooled while the K flux barely 
changed. This can also be seen in DDWS's light curves: af- 
ter JD 2443954, the L and M fluxes grew while K and H 
fell slightly. (DDWS consider their J magnitudes might be 
affected by emission-line flares.) We interpret this phase as 
the continued growth of the newly formed grains while the 
nucleation process itself had slowed; if the grains did not 
grow, then the flux at all the wavelengths would start falling 
at the same time as the grains were carried away from the 
star and cooled as the stellar radiation heating them was 
progressively diluted. After 1979 (DDWS's photometry on 
JD 2444215), fluxes at all wavelengths faded and the colour 
became redder owing to cooling of the dust as the conden- 



sation rate slowed. The CMD shows continued fading, with 
(H — L) reaching a maximum at 1987.2. The fall in aver- 
age dust temperature implied by the reddening in [H — L) 
is interpreted as indicating that the dust being continually 
carried away from the star was no longer being replenished 
by the condensation of new grains at the same rate as be- 
fore. After 1987, the [H — L) gets bluer again while the K 
flux continues to fade. We interpret this as an indication 
that the dust formed in the 1979 maximum had now cooled 
sufficiently that its SED was no longer dominating the L 
band emission. At longer wavelengths, continued cooling af- 
ter 1987 can be seen from the the mid-IR colours (Table [2]): 
e.g. (L' — A^l), which had increased from 1.6 in 1982 to 2.4 
in 1988 (JD 2447519), continued increasing to 2.6 in 1991 
(JD 2448345,7). Some of the SEDs are compared in Fig.[7] 
below. 

Sharp blueward movement in [H — L) accompanied by 
rises in K were observed in 1990 and 1994, the times of the 
mini episodes seen in the light curves. These changes are in- 
terpreted as the consequence of additional episodes of dust 
formation: the bluer {H — L) being a consequence of the 
addition of hot dust at the condensation temperature, while 
the total NIR emission rises. We have sufficiently frequent 
observations to track the evolution of the dust emission in 
the 1994-95 event. This is shown in more detail in Fig.[Ul Af- 
ter the minimum in {H — L), the emission gets redder while 
K continues rising (1994.4-1994.6). This implies that nucle- 
ation of grains has ceased but that the mass of dust contin- 
ues to rise as the grains grow by accretion. After 1994.6, the 
emission fades and the dust continues cooling as it is carried 
away in the wind. Part of a similar clockwise trajectory can 
be seen in 2000.0-2000.6, and smaller bits of the others can 
be discerned. 

After the K minimum in 1997.4, and mini eruptions, 
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Figure 5. Trajectory of WR48a in tlie {H — L), K colour- 
magnitude diagram from the first observation by DDWS in 1978 
to 2011.5. The data are gridded to 0.2-yr bins and the numbers of 
observations contributing to each point are marked with diff'erent 
symbols. 

the {H ~ L) is increasingly influenced by newly formed hot 
dust and remains blue as the K brightens, with a superim- 
posed secondary outburst in 2010.4-2011.3. Generally, how- 
ever, the K, [H — L) slope in the CMD of the major rise to 
maximum differs from those of the mini episodes in that the 
brightening K due to increased dust is not accompanied by 
further blueing of {H — L): the blue limit must be set by the 
dust condensation temperature. 



3.3 Spectral Energy Distributions (SEDs) 

A set of SEDs illustrating the cooling and fading of the dust 
between maximum and our photometry in early 1990 (be- 
fore the mini dust-formation episode described above) de- 
rived from the photometry is given in Fig.[7l Following the 
discussion in Section 1, we de-reddened the photometry with 
a standard interstellar law for Ay ~ 9.0. The dashed line 
through the fluxes at maximum is a 1210-K Planck func- 
tion for reference - but it may be a good description of the 
dust SED at this phase if the emission is optically thick. 
The IR luminosity at maximum (5.6 X 10-"Wm-^ about 
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Figure 6. Detailed view of WR48a in the {H — L), K colour- 
magnitude diagram in 1993—2004 showing trajectories around the 
mini-outbursts in 1994-5, 1999-2000 and, possibly, 1997-8. The 
symbols are as in Fig. [5] and are joined by solid lines when the 
time bins are consecutive and by broken lines when the gaps are 
greater than 0.2 yr. 
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Figure 7. SEDs of WR48a as it faded after maximum in 1979: 
A photometry by DDWS in 1979.5; ffi our mid-IR photometry 
on JD 2445150 combined with DDWS's NIR photometry of 33 
days later; • our photometry in mid-1985; and o our photometry 
in 1990.3. The optical data come from B and V observed in 2006 
by Baume et al. (2009), see text. All the photometry was de- 
reddened by Ay = 9. The dashed line through the 1979 points is 
a 1210-K Planck function for reference. 
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2.8 X 10^1/0 at the adopted distance of 4 kpc) is about 70 per 
cent of the total steUar luminosity, ~ 4 x 10^ for the WC8 
star plus its proposed O companion. This is the fraction of 
the stellar UV-optical flux that is absorbed and re-radiated 
by the dust, and would give the UV-optical depth in the 
dust if it were uniformly distributed around the stars; but 
we know from the images that the dust is not symmeterical, 
so some sightlines will have greater optical depth. 

A Planckian SED in the log XF\ vs. log A plane is iso- 
morphic with temperature and flux level, so the broader 
SEDs seen as the dust cools in 1979-90 indicate that we are 
observing dust having a greater range of temperature. 

We can extrapolate the 1979.5 SED shortwards to 1 mi- 
cron to investigate the possible dilution of the 10830-A and 
9710-A emission lines in DDWS's 1980 spectrum by dust 
emission. This extrapolation, reddened by Ay ~ 9, gives a 
dust-emission magnitude at 1 micron of 10.4. This is sig- 
nificantly brighter than the 1-micron magnitude (~ 12.1) 
extrapolated from the 2006 BVIc photometry by Baume 
et al. (2009) and would provide the additional wavelength- 
dependent continuum emission implied by the dilution of the 
longest wavelength lines in the 1980 spectrum. At shorter 
wavelengths, this is not a problem: the reddened dust- 
emission Ic is 14.8, significantly fainter than the Ic — 13.25 
observed by Baume et al. 

Even at minimum, we can see from the spacing of 
light curves that there was still significant dust emission by 
WR48a. We examine the SED in 1996, near IR minimum by 
combining the SAAO JHKL photometry observed on JD 
2450127, contemporaneously with the /SO-SWS spectrum 
presented in Paper I. Convolution of the SWS spectrum 
with the L filter passband recovers the observed L mag- 
nitude. The SED (Fig. [8]) showing strong dust emission, and 
the absence of further fading after 1997 to the stellar wind 
level, demonstrate continuing dust formation, which replen- 
ishes that carried away in the stellar wind. In this respect, 
WR48a differs from the episodic dust makers like WR140, 
and resembles the variable dust-maker WR 98a discussed in 
Paper I. 

To estimate the SED during the rising branch of the 
light curve, we retrieved mid-IR fiuxes for WR48a observed 
(as AKARI-IRC-Vl J1312398-624255) in the AKARI IRC 
Survey (Ishihara et al. 2009) from the ISAS/JAXA Archive. 
These data come from six scans made between 2006 May 
and 2007 August and we combined them with JHKL ob- 
served in 2006 (JD 2453906 and 2453935). Comparison of 
the AKARI 9-/im and IS-^m fluxes with the JSO-SWS SED 
(Fig.[8| suggests that the mid-IR flux from WR48a barely 
changed over the decade between the two missions, although 
we note that the breadth of the S%W and L18W filters 
used for the IRC survey and their coverage of the 'sili- 
cate' bands requires care in the conversion of in-band fluxes 
to monochromatic fluxes for a SED and the observed and 
extinction-corrected KO-SWS spectra need to be convolved 
with the AKARIYRG optics to examine this further. For the 
present, we note that while the mid-IR flux barely changed, 
the near-IR has increased significantly between 1996 and 
2006, indicating the formation of new hot dust. 

The optical data in Figs [7| and [8] [B, V and Ic) come 
from the photometry in 2006 by Baume et al. (2009). Their 
V = 17.13 agrees weU with the HST WFPC2 magnitude 
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Figure 8. SEDs of WR48a at three phases: A the IR maximum 
1979.5 observed by DDWS; (B close to IR minimum from our 
JHKL photometry in mid-1996 soon after the /50-SWS spec- 
trum from Paper I, solid line; and o while the IR flux was rising 
from our JHKL photometry in mid-2006 with 9-/im and 18-^m 
fluxes observed by AKARI in 2006-07. The optical flux points 
come from BVIc photometry by Baume et al. in 2006, with B 
and V supported by HST WFPC2 43911/ and 556iy photometry 
in 1996, see text. All the photometry was de-reddened by Ay = 9; 
the ISO spectrum was der-reddened by nulling the 9.7-^m silicate 
feature. 

F555VK = 17.08 observed on 1996 July 43 by WaUace et al. 
(2003), while the difference between their B — 19.37 and 
Wallace et al.'s FASdW = 19.79 is readily understood in 
terms of the longer wavelength cutoff of the B filter and 
the red colours of WR48a, which lie outside the range of 
colours for which the transformations between WFPC2 and 
UBV colours given by Holtzman et al. (1995) apply. Hence, 
the B and V are a good measure of the optical flux from 
WR48a at the times of our SEDs in both 2006 and 1996 
at least. We note that the Ic — 13.25 measured by Baume 
et al. is significantly fainter than the i — 12.50 and 12.58 
observed in the DENIS survey (Epchtein et al. 1999) in 1997 
and 1998 respectively. This cannot be emission by the hot 
dust as we can see (Fig.[Sl) that there was more hot dust 
emission in 2006 than in 1997, and further observations in 
this wavelength region are needed to look for variability. 

3.4 Modelling the SED of the 'mini' eruption 

At its simplest, the IR SED of dust emission is given by 
the convolution of the grain emissivity, n\, and the Planck 
function for the dust temperature, Tg-. 

And'^Fx = miHixB{\,Tg) 

where d is our distance to the dust, rrid is the dust mass and 
the emission is optically thin in the IR. Given kx, the tem- 
perature of the grains can be found by fitting the SED over 
a suitable wavelength range, followed by the dust mass from 
the overall fit and distance, d. We calculated kx from the 
optical properties of the 'ACAR' amorphous carbon grains 
prepared in an inert atmosphere in the laboratory by Colan- 
geU et al. (1995) and tabulated by Zubko et al. (1996). Model 

^ Date is from the HST archive. 
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SEDs calculated using the 'AGAR' data gave a better fit 
to broad features in the 5.5-6.5 \x.m ISO spectra of dust- 
making WC8-10 stars (Williams, van der Hucht & Morris, 
1998) than those based on the 'ACH2' grains produced by 
Colangeh et al. in a hydrogen atmosphere - as one might ex- 
pect for grains formed in hydrogen-poor WC stellar winds. A 
comparison of several laboratory analogues of cosmic amor- 
phous carbon grains by Andersen, Lodl & Hofner (1999) 
shows that their absorption efficiencies have fairly similar 
wavelength dependencies, k oc A^^'^, in the IR, but abso- 
lute values ranging by a factor of ~ 5. Consequently, the 
grain temperatures derived from fitting IR SEDs are rela- 
tively insensitive to the choice of grain sample, but not the 
masses of the dust clouds. Usually, circumstellar grains are 
not all at the same temperature, but have a range of temper- 
ature depending on their distances, r*, from the star(s) heat- 
ing them, and also on their size. For small 'AGAR' grains 
and in the absence of significant UV-visual optical depth, 
the temperatures of grains in thermal equilibrium with the 
stellar UV-visual radiation field falls off with distance as 
Tg oc rr°-^* (WiUiams et al. 2009). This dependence is also 
relatively insensitive to grain sample, and was found to ap- 
ply on average even in the full radiative transfer modelling of 
the thick dust spiral of the classical persistent dust-making 
WC9 star WR 104 by Harries et al. (2004). 

Similar techniques will be required for modelling of the 
dust made by WR48a owing to its optical depth. This will 
require knowledge of the dust geometry, and is beyond the 
scope of the present paper; for the present, we consider sim- 
ple optically thin, isothermal modelling of the 1994-5 'mini' 
eruption. We model the SEDs at two epochs: the maximum 
in J and H at 1994.45 (data from JD 2449519) and that 
in K and L in 1994.64 (averaging data of JD 2449583 and 
2449590). This emission is in excess of the continuous 'base- 
line' dust emission by WR48a at the time, and we formed a 
'baseline' SED interpolating between the light curves before 
and after the 'mini' dust-formation episode. The additional 
IR luminosity in the 1994 sub-peak at maximum (i.e. sub- 
tracting the baseline emission) was ~ 10"" Wm-^ equiv- 
alent to ~ 5 X lO'^Z/0 at the adopted distance of 4 kpc, is 
only one per cent of the total stellar luminosity, so we treat 
its emission as being optically thin - although it may suffer 
extinction in the cooler dust in our sightline. Also, if the 
dust is newly formed, the spread in r, and grain size may 
be small enough for the dust to be approximately isother- 
mal, and we can usefully derive an average temperature (Tg) 
from the observed SED. 

The results of fitting these differential SEDs are given 
in Table [71 in which the uncertainties were estimated from 
those in the photometry (0.02 mag. in JHK and 0.03 mag. 
in L). The two epochs roughly separate three phases in the 
evolution of the dust: (a) nucleation of new grains near the 
stars until 1994.45, (b) growth of the grains as they moved 
away from the stars causing the dust mass to increase be- 
tween then and 1994.64 and (c) the continued fading of the 
dust emission. The grain temperature from the 1994.45 fit 
((Tg) = 1220±90 K) is close to that at which the grains 
nucleate; the relatively large uncertainty comes from fitting 
the small difference between the two SEDs, and there must 
be some range in temperature as the first grains to form will 
have moved further from the stars. There was no emission 
from the 'mini' episode in the 1994.06 photometry, so we de- 



Table 7. Isothermal model fits to JHKL SEDs during the 1994 
mini dust-formation episode. 



Date 






formation rate 




(K) 




(Mq yr-i) 


1994.45 


1220±90 


1.4 ±0.7 X 10"* 


3.6 ± 1.8 X 10-8 


1994.64 


1050±20 


4.1 ±0.8 X 10-* 


1.4 ±0.3 X lO-'^ 



termine an average dust formation rate, presumably a com- 
bination of nucleation and growth of the oldest grains, for 
this phase. The lower (Tg) two months later results primar- 
ily from the movement of the mini dust cloud away from the 
stars. There may be a component from more efficient cooling 
by the largest grains, but we do not have the data (e.g. ve- 
locity of the dust grains) to separate the contributions. The 
formation rate in this second phase is about four times that 
in the first and must be almost entirely due to the growth 
of the existing grains as the J and H fluxes were falling at 
the time, implying that nucleation had ceased. The three- 
fold increase in dust mass implies a 40% increase in average 
grain radius in the two months of the second phase. 

The average dust-formation rate, 1.4 ± 0.3 x 10'"^ Mq 
yr-^, seems high for a minor dust-formation episode com- 
pared with the likely mass-loss rate of WR48a (~ 10~^ Mq 
yr-^) and abundance of carbon in its wind, but we recall 
the range in the scales of n\ for different laboratory grain 
samples resulting in a similar range of derived dust mass. 



4 DISCUSSION 

The long-term infrared photometry of WR48a shows com- 
plex behaviour different from that of any other WR star: a 
slow, large-amplitude variation on which are superimposed 
briefer variations of lower amplitude. In all cases in which 
they can be determined, the trajectories in the {H — L), K 
colour-magnitude diagram are loops (clockwise in Figs[S]and 
[Sjl indicative of the formation of new, hot dust, followed by 
its cooling. The recurrence time-scale for the slow variation 
is 32+ yrs; although the 2010-11 NIR photometry reached 
the levels of the earliest (1978-79) pre-maximum photome- 
try by DDWS, the interval cannot be translated into a pe- 
riod owing to the superimposed short-term variations, and 
further observations capturing at least the maximum and 
early decline are needed for a periocJl. Even without a pe- 
riod, there is enough evidence to consider WR48a as a CWB 
periodic dust-maker and discuss it in that context. 

The total dust luminosity at maximum (1979) is 13 x 
that in 1996, near minimum, implying a variable rate of dust 
formation assuming a constant wind dispersing the dust. 
The rate of dust formation is better estimated from the 
short wavelength flux sensitive to the hottest, newly formed 
dust and the range in H between maximum and minimum 
(adjusted for the stellar contribution) gives a factor ~ 16. 
The fraction of the WC8 wind going into the WCR is de- 
termined by the angular size of the WCR, which depends 
on the wind-momentum ratio rj — {A'Ivoo)ob /{Mvoo)wc , 

^ Strictly, observations over another whole cycle, but these would 
take the first author to his centenary and are not planned 
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which does not change around the orbit if the winds collide 
at their terminal velocities, but the pre-shock wind density 
at the WCR does change with the stellar separation if the 
orbit is elliptical. If maximum and minimum dust formation 
coincide with periastron and apastron passages in an ellip- 
tical orbit, and if the dust formation rate is proportional to 
the pre-shock wind density at the WCR, the stellar sepa- 
ration varies by a factor of ~ 4 and the orbital eccentricity 
would be e ~ 0.6, with periastron passages in 1979 and some 
time after mid-2011. 

The mini eruptions appear not to be periodic, arguing 
against the suggestion in Paper I of a third stellar compo- 
nent in WR48a, and other large-scale overdensity structures 
which are intercepted by the WCR as it moves through the 
winds are implied. The only other system having mini out- 
bursts superimposed on its long-term dust-formation varia- 
tion is WR 137. This is one of the few WR stars to have sig- 
nificant continuum polarization and wind asymmetry (Har- 
ries, Hiller & Howarth 1998). Harries, Babler & Fox (2000) 
pointed out the similarity of WR 137 to WR6 and WR 134, 
which display periodic variability in their emission-line mor- 
phology possibly caused by the rotation of large-scale az- 
imuthal structures within their flattened winds. Polariza- 
tion observations of WR48a are needed to to see if this is 
the case here. Also, the Hq emission associated with the 
companion suggests that this might have a flattened wind 
and large-scale structures. 

The persistence of dust formation even at IR minimum 
classifies WR48a as a variable dust maker like WR98a, 
which has a period near 564 d. determined from both the 
rotation of its dust pinwheel (Monnier et al. 1999) and NIR 
photometry (Paper I). The recurrence time-scale of WR48a 
much longer than this, and the periods of the episodic dust 
makers, of which the longest is 13.05 yrs for WR 137, both 
from IR photometry (Williams et al. 2001) and the RV or- 
bit (Lefevre et al. 2005), and even the orbital period, 24.8 
yr, of the persistent dust-maker WR112 inferred from its 
multi-arc dust cloud (Marchenko et al. 2002) . The relevance 
of the long period is the implied size of the orbit (for stellar 
systems of similar mass), with implications for the wind- 
collision process, particularly whether the post-shock WC 
wind is radiative or adiabatic. 

Radiative cooling in the post-shock WC wind is im- 
portant for dust formation (Usov 1991). It was discussed by 
Stevens, Blondin & Pollock (1992), who introduced the cool- 
ing parameter, x'- the ratio of the time-scales for cooling of 
the shocked gas and its escape from the apex of the WCR 
to flow downstream, and given by 

^ = ITr — 

M-7 

where vs is the wind velocity in units of 1000 km s~^, di2 is 
the separation of the stars in units of 10^^ cm, and M_7 is 
the star's mass-loss rate in units of 10^^ Mq yr~^. Stevens 
et al. considered a wind with x ^ 1 to be adiabatic and 
X 1 to be radiative; recent hydrodynamical simulations 
by Parkin & Pittard (2008) find a higher threshold, with the 
compressed wind cooling rapidly if x ^ 3. The WCR in a 
system having a very elliptical orbit, like WR140, can vary 
between adiabatic and radiative as the stellar separation 
varies, becoming radiative close to periastron - the only time 
that dust forms in that system. In the case of WR48a we 



know none of vg, d\2 and M_7, but can estimate x very 
roughly. If the stars in WR48a have masses comparable to 
those in the 7 Vel system they would have d\2 — 530 in a 
P = 33 yr period circular orbit. Adopting a wind velocity of 
1200 km s"^ (Section [23]) and assuming a mass-loss rate of 
lO"'' Mq yr-\ we have x =^ H for the WR48a WCR. This is 
significantly above the threshold and indicates an adiabatic 
WCR - which does not favour dust formation. In an elliptical 
orbit, the separation and hence x would be smaller near 
periastron (e.g. x — 4 for e = 0.6); but they would be larger 
at apastron - and we know from the observations that dust 
forms throughout the orbit. The dust formation by WR112 
raises a similar problem: for the same stellar masses, the 
stars would have d\2 — 390 — 490 in an orbit having e = 
0.11 (Marchenko et al. 2002) and x - 8 - 10 for the same 
estimates of mass loss and velocity, which are also those 
used in their modelling by Marchenko et al. Higher mass- 
loss rates would result in lower values of x, as would lower 
wind velocities, which enter with the fourth power. Wind 
velocities below ~ 900 km s~^ would give values of x — 3 
and possibly radiatively cooling winds. The velocities of both 
stars need to be measured from high-resolution spectra to 
advance the question of the cooling of the post-shock winds 
and the conditions allowing their dust formation. Velocities 
are also needed to aid interpretation of the dust images of 
the two stars. 

The recurrence of dust formation, the dilution of the 
WC emission lines implying presence of a luminous com- 
panion, the anomalous emission in Ha (and possible H/3 
and Br7) supporting this, and the strong X-ray emission re- 
ported by Zhekov et al. (2011) and Mauerhan et al. (2011) 
all point to WR48a being a long-period CWB. The IR pho- 
tometry needs to be extended to determine a period and 
characterise the 'mini' outburst and the spectrum needs to 
be re-observed at higher resolution and S/N to elucidate the 
types of both components and determine the WC wind ve- 
locity to further the consideration of the confiict between 
observed dust formation and an apparently adiabatic post- 
shock wind. Knowledge of the wind cooling will also be valu- 
able in understanding the X-ray flux from WR48a and its 
long term variation. The long period of the system suggests 
that any non-thermal radio emission from the WCR will 
have little difficulty in escaping the wind (cf. Dougherty 
& Williams 2000). Its distance is comparable to that of 
WR112 (4 kpc, van der Hucht 2001), which shows vari- 
able non-thermal radio emission at the level of ~ a few mjy 
(Chapman et al. 1999, Monnier et al. 2002), so we expect 
the similar fiux levels from WR48a if the stars have simi- 
lar winds, but with greater amplitude of variation given the 
greater eccentricity and inclination of the WR48a system. 
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